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We study the effect of antiferromagnetic interactions on the single spin-flip Glauber dynamics of 
two different one-dimensional (ID) Ising models with spin ±1. The first model is an Ising chain 
with antiferromagnetic exchange interaction limited to nearest neighbors and subject to an oscil- 
lating magnetic field. The system of master equations describing the time evolution of sublattice 
magnetizations can easily be solved within a linear field approximation and a long time limit. Res- 
onant behavior of the magnetization as a function of temperature (stochastic resonance) is found, 
at low frequency, only when spins on opposite sublattices are uncompensated owing to different 
gyromagnetic factors (i.e., in the presence of a ferrimagnetic short range order). The second model 
is the axial next-nearest neighbor Ising (ANNNI) chain, where an antiferromagnetic exchange be- 
tween next-nearest neighbors (nnn) is assumed to compete with a nearest-neighbor (nn) exchange 
interaction of either sign. The long time response of the model to a weak, oscillating magnetic field 
is investigated in the framework of a decoupling approximation for three-spin correlation functions, 
which is required to close the system of master equations. The calculation, within such an approx- 
imate theoretical scheme, of the dynamic critical exponent z, defined as 1/r ~ (l/£) z (where r is 
the longest relaxation time and £ is the correlation length of the chain), suggests that the T = 
single spin-flip Glauber dynamics of the ANNNI chain is in a different universality class than that 
of the unfrustrated Ising chain. 

PACS numbers: 75.10.-b, 75.10.Pq, 75.50.Ee, 75.50.Gg 

I. INTRODUCTION 

After the publication of fundamental papers!*- 2 on stochastic resonance (SR), it was realized that the response 
amplitude of a nonlinear dynamic system to an external periodic signal is greatly enhanced as a function of noise 
strength, in the presence of a matching between the frequency of the external force and the escape rate across an 
intrinsic energy barrier. Most of the SR research^ was pursued on dynamic systems with a double well potential, 
subject to both periodic and random forces, while only a few investigations of SR in extended or coupled systems 
have yet been conducted^ 

The Ising model with Glauber dynamics 5 - can be viewed as a set of coupled two-state oscillators, where the coherent 
signal is provided by an external oscillating magnetic field and thermal fluctuations are the only source of random 
noise. Each spin is assumed to be in interaction with a heat reservoir of some sort, which causes it to flip between 
the values a = +1 and a = —1 randomly with time. In the presence of magnetic coupling between the spins, the 
transition probability for one spin to flip is assumed to depend on the configuration of the neighboring spins. The 
time evolution of the system is described by a master equation where the transition rates verify the detailed-balance 
condition. Solving the master equation, the time dependence of the magnetization and of the spin correlation functions 
can be obtained. For exchange interaction limited to nearest neighbor (nn) spins, the response of the Ising model with 
Glauber dynamics to an oscillating magnetic field was investigated in one (1D)^^ two (2D)) 7 - and three (3D)&2 spatial 
dimensions. For the ID nn Ising ferromagnet, Brey and Prados- obtained an analytic expression, within the linear 
field approximation, for the amplitude and the phase of the induced magnetization. The amplitude always presents 
a maximum as a function of temperature, with a genuine resonant behavior only for low frequencies. The Glauber 
dynamics of the ID Ising model with antiferromagnetic next-nearest neighbor (nnn) exchange interaction competing 
with the nn one was investigated by Yang, 10 - who employed a decoupling approximation to solve the master equation 
and get an analytical expression for the time-dependent magnetization. He also found, by heuristic arguments, the 
dynamic critical exponent z, defined as 1/r « (I/O 2 (where r is the longest relaxation time and £ is the correlation 
length of the chain)ii, to be z — 2, the same as that of the unfrustrated ID nn Ising model. 

In this paper, we study - at finite temperature T > - the effect of antiferromagnetic (AF) exchange interactions on 
the single spin-flip Glauber dynamics of two different one-dimensional Ising models. Our interest in kinetic ID Ising 
models with AF interactions is motivated by recently sinthesized cobalt-based 12 J=S and rare-earth-base d 14 ' 15 single 
chain magnets, showing slow relaxation of the magnetization at low temperature. The magnetic properties of the 
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former chain compound, [Co(hfac)2NITPhOMe], can be described in terms of a ID Ising model with AF nn exchange 
coupling. 13,16 However, the resulting short range order is ferrimagnetic, owing to the alternation along the chain of two 
different kinds of magnetic centers (a metal ion, Co 2+ , and a nitronyl-nitroxide radical, PhOMe), both with S = 1/2 
but with different gyromagnetic factors. In spite of further complications due to non-collinearity of the spins ji 6 - this 
system was shown to be the first experimental realization of a ID nn Ising model with Glauber dynamics.-! 3 " The single 
chain magnets belonging to the latter class of rare-earth-based compounds, of general formula [M(hfac)3(NiTPhOPh)], 
where M=Eu, Gd, Tb, Dy, Ho, Er, or Yb, and PhOPh is a nitronyl-nitroxide radical, are characterized by strong Ising- 
type anisotropy and by the simultaneous presence of both nn and nnn exchange interactions between the magnetic 
centers, with the last ones being antiferromagnetic in naturei 14 i 15 

The paper is organized as follows. In Section II we investigate the Glauber dynamics in a collinear Ising chain, 
with antiferromagnetic exchange interaction limited to nearest neighbors and different gyromagnetic factors on the 
two opposite sublattices, subject to an oscillating magnetic field. The system of master equations describing the 
time evolution of sublattice magnetizations can easily be solved within a linear field approximation and a long time 
limit. Resonant behavior of the magnetization as a function of temperature (stochastic resonance) is found, at 
low frequency, only when spins on opposite sublattices are uncompensated owing to different gyromagnetic factors 
(i.e., in the presence of a ferrimagnetic short range order). In Section III we investigate the ID axial next-nearest 
neighbor Ising (ANNNI) model, where an antiferromagnetic exchange between next-nearest neighbor spins is assumed 
to compete with a nearest-neighbor exchange interaction of either sign. The long time response of the model to a 
weak, oscillating magnetic field is investigated in the framework of a decoupling approximation (required in order to 
close the system of master equations) for three-spin correlation functions, which in principle is more accurate than the 
one reported in Ref. ITol . As a consequence, our approximate calculation of the dynamic critical exponent z suggests 
that the T = single spin-flip Glauber dynamics of the ANNNI chain is in a different universality class than that of 
the unfrustrated Ising chain. Finally the conclusions are drawn in Section IV. 



II. GLAUBER DYNAMICS IN THE NEAREST-NEIGHBOR FERRIMAGNETIC ISING CHAIN 

We consider a one-dimensional Ising model with a nearest-neighbor antiferromagnetic exchange interaction, J < 0, 
in the presence of a time-dependent external field. The Hamiltonian of the system is 

N N/2 

U = - JJ2 VjVj+i - Mo H(t) ]T(<m ol 3 _ x + g B a* 2j ) (1) 
i=i i=i 

where /zo is the Bohr magneton, and H (t) = Hoe~ luJt is an external magnetic field applied along the z direction and 
oscillating in time with frequency uj. Spins on opposite sublattices are allowed to take possibly different gyromagnetic 
factors (gA ^ 9b), while we assume erj = ±1 Vj. Hereafter, the z index will be dropped for ease of notation. In the 
absence of a magnetic field, if gA ^ 9b the ground state is ferrimagnetic, with opposite uncompensated magnetizations 
on the two sublattices; if gA — gB the ground state is antiferromagnetic, with compensated sublattice magnetizations. 

When the system is endowed with single spin-flip Glauber dynamics^ its time evolution is described by the master 
equation 

^p(a, t) = [WiiRjtrMRj*, t) - W s (a)p(a, f)] (2) 
j 

where p(a, t) is the probability for the system to assume the configuration a = {ci, • • • , aj, • • • , <jn} at time t, Rj<J 
is the configuration obtained from a by flipping spin j, and Wj(a), Wj(Rja) are the transition rates between such 
configurations. 

For a ID Ising model of spins [(Xj = ±1) with ferromagnetic nn exchange interaction J > and gyromagnetic factor 
g, Brey and Prados 6 showed that, for low frequency, a stochastic resonance phenomenon occurs: i.e., the induced 
magnetization M(t) — gpo^2j=i(crj',t) oscillates at the same frequency as the magnetic field, and the amplitude of 
M(t) presents a sharp maximum as a function of temperature T. The resonance temperature, T T , is determined by 
the matching between the frequency to, of the external field and the inverse of the statistical time scale, l/r(T r ), 
associated to the spontaneous (i.e., in zero field) decay of the magnetization. In zero field, the magnetization of the 
ID nn Ising ferromagnet was found^i 6 - to relax to its equilibrium value, M = 0, with the asymptotic t — > oo behavior 
M(t) ~ e~*/ T ( T '/\/L The relaxation time r was found to be exponentially divergent for T — > 0, r(T) w e^r (where 
k denotes Boltzmann's constant), and to become of the order of the inverse of the transition rate of an isolated spin 
forT^oo, r^l/a^. 
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For a ID Ising model with antiferromagnetic nn exchange interaction (J < 0), the master equation |(3J) is still the 
starting point for the study of the chain dynamics. In this case, if g A ^ <7b, the transition rates in the presence of a 
field are assumed to be different for even (A) and odd (B) lattice sites j 



Wj{a) = Wf\a) 1 - aj tanh(^,B) 



1 

2 a 



1 - <7j tanh(/3 a,b) 



(3) 



where W- (cr) denote the transition rates in zero field. The transition rate of an isolated spin, |a, is considered as 
temperature independent and sets the time scale. In the case of interacting spins, the probability per unit time of the 
j-th spin to flip depends on the orientation of its nearest neighbors. The magnetic field favors one orientation with 
respect to the other. A correspondence between the parameters 7, Pa,b of the stochastic model and the parameters 



J, gA,B^oH{t) of the statistical Ising model can be obtained^ observing that at equilibrium ^p(a,t) 
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Next, requiring the detailed balance (i.e., the microscopic reversibility) condition to be satisfied 

WjiRjo) _ Peq {<r,t) 



Wj(a) p eg {Rj(T,ty 



with p eq (a,t) 



e fcT 



and p eq (Rj(T, t) 



e kf , one readily obtains 



tanh(^), 



Am 



tanh 



9A,B^oH(t)\ 



kT 



0, so that 
(4) 

(5) 

(G) 



The evolution equation for the spin expectation value (&j;t) ~ <jjp(a,t) is directly obtained from the mas- 
ter equation to be j^(<7j;t) = —2(ajWj(a);t)}^ Considering that for model ([T]) the spins belong to two opposite 
sublattices, the system of evolution equations in the presence of an oscillating field is 

d 



d(at) 
d 



d(at) 



(a 2j -i;t) = -(a 2] -i;t) + ^j^(<r 2 j-2;t) + {a 2j ;t)^ + (3 A 1 - ^7 ((^-2^-1; t) + (a 2] -i<7 2] ; i)) 
(<j 2j ;t) = -(a 2j ;t) + ]p({<J 2] -i\t) + (a 2j+1 ;t)j + f3 B 1 - ^l({cr 2 j-ia 2j ;t) + (a 2j a 2 j+i;t)^ 



(7) 



The system is not closed owing to the presence of two-spin, time-dependent correlation functions on the right 
hand sides. In order to solve it, a linear field approximation is made^ so that taj3h[(gA, B ^oHo)/(kT)] can be 
expanded for small values of the argument and two-spin correlations can be evaluated in the absence of a field. 
Moreover, if in the long time limit t — > 00 the nn correlation functions are assumed^ to take their equilibrium value 
(cTjCj+i; t) — ► T] = tanh[J/ (kT)], the system of two coupled equations of motion for the two sublattice magnetizations 



M 1 {t)=g A iM) 



N/2 

E 



(<j 2] -i;t), M 2 (t) 



N/2 
3=1 



(8) 



can be written in matrix form 

' 9 Mt) 



d(at)- 
d 

mat) 



M 2 (t) 



1 -24-7 

9b 1 

9A ' 



(9) 



Taking into account that 7 = rrr? , the temperature dependent coefficient j\f(T) can be expressed as 



A^) = ^^(l-7.) = -^^4. (10) 

V ' 2 kT y ' U 2 kT 1 + rf V ' 

The above system can be decoupled diagonalizing the 2x2 non-symmetric matrix on the l.h.s. of Eq. Denoting 
by A4i(t) and M 2 (t) the normal modes, one obtains 



M 2 (t) r \ X 2 J V M 2 (t) J J *^'\f 2 
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where the eigenvalues A rl (n = 1, 2) turn out to be independent of the gyromagnetic factors gA and gs 

Ai = 1 - 7, A 2 = 1 + 7, 

and the /„ (n = 1, 2) coefficients are 



/i 



(.9b+.9a) /2 



,9s 



(.9s -9a)- 



2 w " y " 2 

The relationships between the normal modes A4 n (t) and the sublattice magnetizations M n (t) (n — 1, 2) are 



Mi(t) 



M 2 (t) 



5s 

9A 



Mi(t) 



M a (t) = - Af a (t) 



SB 

.9/1 



Mi(*) 



(11) 



(12) 



(13) 



(i.e., A4i(i) and M 2 (t) are related to the net and the staggered magnetization, respectively). Conversely, one has 

9A 



Mi(t) = Mi(t)-M 2 (t) 
9b L 

The general solution for the normal modes is (n = 1, 2) 



M 2 {t) = Mi(t) + M 2 (t). 



M n {t)=M n (to)e-— +N{T)f n dt'e 



(14) 



(15) 



Co 



where the relaxation times t„ are expressed, in terms of the eigenvalues A„ of the non-symmetric 2x2 matrix, as 
t„ = l/(a\ n ), so that 



1 



1 



n 



T2 



(16) 



a(l — 7) ' a(l + 7) 

In the absence of an external magnetic field, J\f(T) — 0, the normal modes M n (t) are found to relax exponentially. 

In the low temperature limit, T — > 0, one has 7 = tanh (fy) ~ ]7](^ ~ 2e - "^^, so that for antifcrromagnetic 

nn exchange (J < 0), the first relaxation time is simply t\ « ^j, while the second relaxation time is exponentially 

diverging with decreasing T, t 2 ~ 2a e "^' F° r hiS n temperatures, kT 3> both relaxation times become of the 
order of the inverse of the transition rate of an isolated spin, t\ k, t 2 w 1/a. 

For non vanishing magnetic field, the time dependence of the normal modes is obtained letting to — * —00 in Eq. (|15p 



M n (t) = Af(T) 



fn 



A„ 1 



(n = 1,2) 



The total magnetization is 

M tot (t) = M 1 (t) + M 2 (t) 



9B + 9A 
9B 



Mx{t) 



gB - gA 

9B 



M 2 {t) = X (u,T)H e- 



where the complex susceptibility x(u,T) is given by 

'9b + 9A\ 2 l+V 



1 — ry 1 — iuTi 



9b - 9A\ 2 1 - f] 



1 



1 + rj 1 — iuiT 2 



(17) 



(18) 



(19) 



In the limit u> — > 0, the static susceptibility of the Ising ferrimagnetic chain in zero field is correctly recovered: 
see Appendix A.l, Eq. (1A6|) . for details. As regards the dynamic response of the system to a weak, oscillating 
magnetic field, from Eq. ([IT))) it is apparent that, for antiferromagnetic nn exchange (J < 0) and T — > 0, the first 
term on the r.h.s. is associated with a fast relaxation, while the second term with an exponentially slow relaxation. 
Thus, a resonant behavior, similar to the one observed in the ferromagnetic nn Ising chain endowed with single spin- 
flip Glauber dynamics,— is possible only when spins on opposite sublattices are uncompensated owing to different 
gyromagnetic factors (i.e., in the presence of ferrimagnetic short range order). See Fig. 1, where the temperature 
dependence of the amplitude of x(w,T) is reported, for selected values of the frequency, both in the compensated 
( J < and gA = 9b) and uncompensated ( J < and gA ^ gs) case. 

The resonant behavior shown by the ferrimagnetic chain at low frequency (see Fig. lc) is a manifestation of the 
stochastic resonance phenomenon^, i.e., the response of a set of coupled bistable systems to a periodic drive is 
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kT/J kT/J 

FIG. 1: (color online) Temperature dependence of the amplitude of the complex susceptibility |x(w,T)| for an Ising chain 
with antiferromagnetic nearest neighbor interaction J = — 1, subject to a weak external magnetic field oscillating at frequency 
uj. Figures (a), (b) refer to the compensated case {gA = Qb = 2), while Figures (c) and (d) to the uncompensated case 
(gA = 2, gs — 3), for selected values of the frequency (u)/a = 0.001 and 10). In Figs. (c,d) the thin (color) lines represent the 
contributions to the amplitude of the two terms on the r.h.s. of Eq. (|19[) . while the thick (black) line is their sum. A resonant 
behavior (similar to the one predicted for the nn Ising ferromagnetic chain endowed with single spin-flip Glauber dynamics)^ 
is observed only in the uncompensated case for low frequency (notice the enhanced vertical scale in Fig. c). 
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FIG. 2: (color online) Frequency dependence of the peak temperature of the amplitude of the complex susceptibility |x(w,T)| 
of an Ising chain with nearest neighbor exchange interaction. Triangles: compensated antiferromagnet (J = —1, gA — gs = 2); 



circles: uncompensated ferrimagnet (J 
lines are guides to the eye. 



-1, gA — 2, gB — 3); squares: ferromagnet (J = +1, gA = gB = 2).— The dashed 
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enhanced in the presence of a stochastic noise when a matching occurs between the fluctuation induced switching 
rate of the system and the forcing frequency. In the ferrimagnetic chain, the role of stochastic noise is played by 
thermal fluctuations and the resonance peak occurs when the deterministic time scale of the external magnetic field 
matches with the statistical time scale associated to the spontaneous decay of the net magnetization M to t (t) ■ For low 
frequency lu -C a (i.e., low temperature), the resonance condition for the uncompensated case is 



UJ 1 PS T 2 (Tp eak ), 



(20) 



while for the compensated case only the mode with fast relaxation n ps 0{a~ 1 ) contributes, providing a broad peak 
rather than a genuine resonance. For high frequency lu 3> a (i.e., high temperature) a broad peak is found, both for 
the uncompensated and the compensated case, since the two relaxation times T\ and t 2 become of the order of 1/a, 
so that the resonance condition cannot be fulfilled £ 

The frequency dependence of the peak temperature T peak is reported in Fig. 2 both for the compensated (anti- 
ferromagnetic) and the uncompensated (ferrimagnetic) chain, and compared with the ferromagnetic counterpart. 6 
In the compensated case, the frequency dependence of the peak is very smooth, owing to the smooth temperature 
dependence of the relaxation time t%, ranging between l/(2a) at low T and 1/a at high T. In the uncompensated 
case, a behavior very similar to the ferromagnetic one is observed for low frequency: the reason is that for low lu the 
dominant contribution to x( w ) is provided by the second term on the r.h.s. of Eq. (|19p . At intermediate frequency, 
a maximum is observed owing to the coming into play of the first term on the r.h.s. of Eq. (|19[) . Finally, for lu 3> a, 
the amplitude of x(w) becomes 



\ x{l u,T)\*N§° 
kl lu 



(££±-)^ (1 _ 7) + ( £^) 



9b - 9A\ 2 1-r) 



1 



V 



(1 + 7) 



(21) 



where both terms in square brackets on the r.h.s. of Eq. (|21[) present a maximum at the same temperature, which is 
numerically determined to be T pea k ps 1.66711|J|. 



III. GLAUBER DYNAMICS IN THE AXIAL-NEXT-NEAREST-NEIGHBOR-ISING (ANNNI) CHAIN 

We consider a ID axial-next-nearest-neighbor Ising (ANNNI) model with spins alternating on two interlacing 
sublattices (denoted by A and B), with Hamiltonian 



N/2 N/2 N/2 

n = Jl V (<4i-l< + -J^Y, + CT 2>2 1+2 ) - fM)H(t) (gA<4i-l + g B <4,) ■ ( 22 ) 



>1 

1=1 



The intra-sublattice antiferromagnetic next-nearest neighbor coupling J 2 < competes with the inter-sublattice near- 
est neighbor coupling J\, which may be of either sign. In what follows, we shall assume J\ > (ferromagnetic 
coupling). H(t) — Hoe luJt is an external magnetic field applied along the z direction and oscillating in time with fre- 
quency lu, fiQ denotes the Bohr magneton and the spins erf = ±1 are allowed to assume possibly different gyromagnetic 
factors on odd and even sites (gA ^ <?s); the z index shall be dropped for ease of notation. 

In the limiting case <?a = 9b = g, Eq. (1) reduces to the well-known ANNNI (axial next-nearest-neighbor Ising) 
model.^ 7 Depending on the competition ratio r = — J2/J1, this model in zero field is known to admit a ferromagnetic 
ground state for r < 1/2, and a (2,2) antiphase structure (two spins up, two spins down), with zero magnetization, 
for r > 1/2; for r = 1/2 the ground state is degenerate and disordered.— At finite temperatures, the model cannot 
support long range order; however, a strong short range order is present in the paramagnetic phase. For zero applied 
field, as far as the thermodynamic properties are concerned^ the ID ANNNI model can be mapped into an equivalent 
ID Ising model with only nearest neigbor interaction in an effective field, and analytic results (see Appendix A. 2) can 
be obtained for the partition function and the spin correlation functions * 21 ! 22 In the presence of a static magnetic field, 
the ground state of the generalized ANNNI model, i.e. a chain of alternating spins with different quantum numbers 
and different nnn exchange interactions on the two sublattices, was thoroughly investigated? 2 ^ and the thermodynamic 
properties were exactly calculated (though numerically) by the transfer matrix methodj 23 ' 24 

Here we aim at investigating the long-time dynamic response of the ANNNI chain, Eq. (|2"2"|) , to a weak, external 
magnetic field oscillating in time. The time evolution of the system is still described by the master equation @, but 

with respect to the case of the nn Ising chain, the transition rates in zero field, Wf\a), are now assumed to take the 
form 



W,- V) = 7; a 1 - t;7iCj0j-i + 1 - ^72^(^-2 + a j+2 ) (23) 



1 



2 



1 



2 



1 



2 
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meaning that the probability per unit time of the j'-th spin to flip depends on the status of both its nearest neighbors 
and next nearest neighbors; ^a, the transition rate of an isolated spin, is arbitrary and sets the time scale. In the 
presence of a field applied along the z axis, the transition rates Wj(a) are given by 



Wj(a) = Wj 0) (a) \l - aj tanh(/3 A , B ) 



(24) 



As usual, a correspondence between the parameters 71, 72, Pa,b of the stochastic model and the parameters J\, 
J2, 9A,BHoH(i) of the statistical ANNNI model can be obtained requiring the detailed balance (i.e., the microscopic 
reversibility) condition, Eq. ([5]) to be satisfied at equilibrium. One findsiH 



71 = tanh 



-tt^I < 72 = tanh ( — 
kT J ' ' V kT 



11 



A,B 



tanh 



' 9A,B^H(t) 

kT 



(25) 



The stochastic equation of motion for the spin expectation value (o~j;t) = o~jp(o~, t) in the presence of an 
oscillating field is then obtained, from the master equation, to be J^Ojjt) = — 2(<7jWj(a); t), giving^ 



0_ 

Of 



(cr-t) = -(<jf,t) + ^7i((fj-x;*) + {a-j+i;t)j + -72^-2;*) + (o-j+2;*)) 

- ^7l72((c r jCr i _iCr i _ 2 ;i) + {(TjO-j-l<Tj+r,t) + {aj<Tj + l<Tj-2;t) + (<Tjaj+l(Tj + 2;t)^ 

9A,B^oH(t) ^ 



tanh 



kT 



x 1 - ^7i((o-iOj-i;*> + (ojOj+i;*>) - 7p {i,<7 3 <Tj- 2 \t) + (ajcr j+2 ;tf) 
+ ^7i 72 {[(Tj-i(T 3 -2\t) + (aj-i<Tj +2 ;t) + (a j+1 aj- 2 ;t) + (<T j+1 <T j+2 ;t)) 



(26) 



where we remind that the subscripts A and B refer to the case of j odd and j even, respectively. This set of equations 
is not closed, owing to the time-dependent two-spin and three-spin correlation functions on the r.h.s. In order to solve 
it, we make the following approximations. 

• For sufficiently weak fields (x = (gA,B^oHo)/(kT) <§; 1), the hyperbolic tangent on the r.h.s. of Eq. (j26|) is 
expanded for low values of the argument (tanh x w x) and two-spin correlation functions are calculated in the 
absence of a field. 

• Three-spin correlation functions are decoupled, in all possible ways, into products of a single-spin expectation 
value and a two-spin correlation function 



(ajO-j +m o-j +n ;t) w (aj;t){o-j +m aj +n ;t) + (a j +m ;t)(o a j +n ;t) + (aj +n ;t){ajaj +m ;t). 



(27) 



Notice that a different, and incomplete, decoupling was adopted in Ref. [l0|, thus leading to different results with 
respect to the present work. 

• For sufficiently long times, two-spin correlation functions between n-th neighbors are assumed to be independent 
of the initial conditions and to take their static equilibrium values (o~j<jj+ n ;t) — > r\ n for t — > 00. Static two-spin 
correlation functions r] n = (aj<jj+ n ) can be exactly calculated in ID via the transfer matrix methodi 21 i 22 ' 23 i 24 
For qa = 9b and Hq = 0, analytic result a 21 ' 22 can be obtained for rj n : see Appendix A. 2 for details. 

Under these approximations, the master equation for the spin expectation value on a generic site j becomes 



1 + ^7172(^1 + m) (°j;t) 



+ ^li(^--l2m) (<rj-i',t) + (<Tj+i;t) +-72(l-7i?7i) {o-j-2]t) + (o- J+2 ;t) 
\gA,BHoH(t) 



kT 



1 - 71V1 - J2V2 + tj7i 72(?7i + 



(28) 



In the range of the competition ratio r corresponding to weak nnn antiferromagnetism < r < ^ , the ground state 
of the model is ferromagnetic (since we have assumed J\ > 0), while for strong nnn antiferromagnetism | < r < 00, 
the ground state is the so-called (2, 2) antiphase state, consisting of two spins up followed by two spins down. The 
two different regimes shall be investigated separately since they require different order parameters. 
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A. Weak nnn antiferromagnetism (competition ratio < r < |) 



In the range of the competition ratio r corresponding to the ferromagnetic ground state (0 < r < j), owing to the 
different gyromagnetic factors on odd (gA) and even (gs) lattice sites, it is necessary to consider the magnetizations 
over two sublattices, like in Eq. ([5]), as the order parameter. From the master equation, Eq. (|2"8"j). one is thus led to 
consider a system of two coupled equations of motion, which can be written just like Eq. @, with the elements of the 
2x2 non-symmetric matrix now being 



1 



an 

012 



1 - 72(1 - 7i»7i) + 7:7i72(?7i + = 



9A s 

71 (1 - 72?7 2 ) 

9B 



2 

021 



022 



71 (1 - 72772) 

9A 



and the temperature dependent coefficient 

y ' 2 kT 

After diagonalization, the eigenvalues now turn out to be 



1 

1 - 7i7/i - 727/2 + 777271(771 



1 



773) 



(29) 



(30) 



Ai = 1 - 71(1 - 72772) - 72(1 - 71771) + 27172(771 + m) 

A 2 = 1 + 7i(l - 727/2) - 72(1 - 7i77i) + 777172(771 + %), 



(31) 



independent of the gyromagnetic factors gA and gs- The relationships between the normal modes M n (t) and the 
sublattice magnetizations M n (t) (n = 1,2) are the same as in Eqs. (TIB")) , Also the expressions for the /-coefficients 



are the same, i.e. fi — ^f-(gB + gA), f% = ^(gB — gA)- As before, the general solution for the normal modes takes 
the form in Eq. (|15[) . where the relaxation times are r n — 77^—, with the eigenvalues now given by Eq. (|31j) . Finally, 
in the case of weak nnn antiferromagnetic coupling, the complex susceptibility of the ANNNI chain turns out to be 



Mo 
kT 
9b 



(^) 



1 - 7x771 
9A\ 2 1 



727?2 
1 



:7l72(77l +7/3) 



9b - gA\ 2 1 



Ai 1 — iu)T\ 



1 



A2 1 — 1UJT2 



(32) 



In the limiting case r = 0, the well-known result for the nn Ising chain 5 '- is correctly recovered. In the case 
< r < ~, we show in Fig. 3 that the approximate static susceptibility, calculated from Eq. (|32p for zero frequency, 
turns out to be in good agreement with the exact transfer matrix result ) 21 ! 22 Eq. (|A14[) . only at high temperatures 
(kT > Ji). In contrast, an unphysical (negative) static susceptibility is obtained at low temperatures, as a consequence 
of the negative values assumed by the eigenvalue Ai for kT < J\. 

The low-temperature failure of Eq. (|32p can be attributed to the decoupling (f2"T|) of three-spin correlation functions, 
which was made in order to close the set of master equations, Eq. (|26]) : in fact, decoupling approximations have the 
drawback to be uncontrollable, but in principle they are expected to be more accurate the higher the temperature. 
Moreover, at low temperatures one can guess another source of error to lie in the assumption that, for sufficiently long 
times, the spin-spin correlation functions take their static equilibrium values: (ajaj-\- n ; t) — > rj n for t — > 00. In fact, 
for competition ratio in the range < r < 1, the ID ANNNI model with Glauber dynamics is known to be lacking 
in ergodicity at T = 0: the ground state can not be reached by single spin-flip Glauber dynamics, after a sudden 
cooling of the system down to T — starting from high temperature. The difference between the static (r = 1/2) 18 
and the dynamic (r = l) 2 ^ ground state phase boundary of the ID ANNNI model was pointed out by Redner and 
Krapivsky 2 ^, who showed that for < r < 1/2 the ferromagnetic ground state can not be reached because of the 
repulsion between domain walls which forces them to be at least two lattice constants apart, while for 1/2 < r < 1 
the (2, 2) antiphase ground state can not be reached owing to the persistence of isolated domains of length > 33£ In 
contrast, both for r = (ID nn Ising model)S.i2& and r > 1 (ID ANNNI model with strong nnn AF coupling) 2 — the 
ground state can asymptotically (t — > 00) be reached at T = 0. 

The low temperature failure of our approximate theory in the case < r < h prevented us from calculating the 
temperature dependence of the amplitude of the complex susceptibility. However, it is worth observing that, since for 
T — > the zero-field static susceptibility diverges, 21 a resonant behavior might be expected for low frequency provided 
that the system admits also a diverging relaxation time for low temperature. 
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2 4 6 8 

kT/J, 

FIG. 3: (color online) Temperature dependence of the static susceptibility x( w = 0,T) of an ANNNI chain with Ji = 1, 
Ji = —0.35 and gA = Qb = 2, corresponding to a value r — 0.35 of the competition ratio (weak nnn antiferromagnetism) . The 
thick line is the exact transfer matrix result, while open circles denote the approximate calculation, Eq. (|32p . The temperature 
dependence of the eigenvalues Ai and A2 (|31[) is also shown by the dashed lines. The approximations made to close the set of 
master equations (|26[l are found to fail for low temperatures. 



B. Strong nnn antiferromagnetism (competition ratio | < r < 00) 

In the range of the competition ratio r corresponding to the (2, 2)-antiphase state (i < r < 00), it is necessary to 
consider the magnetizations over four sublattices^l 

N__-i 

4 L 4 1 

Mi(t) = #aMo ^2 (ci+4j-;i), M 2 (t) = gsHo ^2 (o- 2 +4j;t) 



y -1 , 

4 1 4 1 

M 3 {t) = gAUo ^2 (a 3+ 4j;t), M 4 (t) = g B ^o (<r 4+4j ;t) 

j=o j=o 



(33) 



as the order parameter. One is thus led to consider a system of four coupled equations of motion, which can be 
written in matrix form as 



where 



and 



-M- 



(t)\ 



d(at 



A 



A B C B\ /Mi(t) 

V A V C M 2 {t) 
C B A B \\ M 3 (t) 

V C V A J \M A (t) 



- N(T) 




1 



an = a 2 2 = a 33 = a 44 = 1 + -7271(^1 + %) 

B = a 12 = a 14 = a 32 = a 34 = -- — 71(1 - 72^2) 

2.gs 

C = ais = a 31 = a 24 = a 42 = -72(1 - 71^1) 

L> = a 2 i = a 23 = 041 = a 43 = -- — 7i(l - 72772) 

2 gA 



Af(T) 



4 kT 



1 - 7i»?i - 72?72 + 27271(^1 + m) 



(34) 
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FIG. 4: (color online) Temperature dependence of the static susceptibility x( w = 0,T) of an ANNNI chain with Ji = 1, 
gA = 9b — 2, for two different values of the nnn exchange constant: (a) J2 = —0.75 and (b) J2 = —1.25, corresponding to 
competition ratio | < r < 1 and r > 1 respectively (strong nnn antiferromagnetism) . The thick line is the exact transfer 
matrix result, while open circles denote the approximate calculation, Eq. (|4ip . The temperature dependence of the eigenvalues 
Ai = A2, A3 and A4 ((36)) is also shown by the dashed lines. 



Diagonalizing the matrix of coefficients, the time dependence of the eigenmodes is found to be (n = 1, 2, 3, 4) 

M n (t) = M n (t )e-^r +Af(T)f n f dt'e^e^' (35) 

to 



where r n = are the relaxation times and fx = fi = 0, fa = ^-{gB — 9a), fi = ^{qb + 9a)- The eigenvalues 
(A„) of the 4x4 nonsymmetric matrix of the coefficients turn out to be independent of §a and 93 



Ai = A 2 = 1 +72(1 - 7i Vi) + 7^7271 {vi + m) 



11 



a 3 = i +71(1 - 72^2) - 72(1 - jim) + ^7271 (^i + m) 
A 4 = 1 - 71(1 -72*72) - 72(1 - 71V1) + 2^271(^1 +m) 

For non vanishing magnetic field, the time dependence of the eigenmodes is 

In I • 



M n (t) = 7V(T)- 



(n= 1,2,3,4) 



A„ 1 - iujT n 

The relationships between the eigenmodes A4 n (t) and the sublattice magnetizations M n (t) (n = 1,2,3,4) are 



(36) 



(37) 



and conversely 



Mx{t) = \ 

M 2 (t) = l - 

M 3 (t) = \ 
1 
4 



Mi(t) 



M 4 (t)-M 2 {t) 
M a (t)-Mi(t) 
M 4 (i)+M 2 (i) 



.9b 



M 3 (t)+Mi(t) 



M 4 (t) + M a (t) + — M 3 (t) + Mi(t) 
9A 



(38) 



Mi(t) = ^ M 4 (i - M 3 *) - -Ma (t) 

5b L J 
M a (t) = Mi(t) + M 3 (t) -Mi(t) 

M 3 (t) = — \M 4 (t) - M 3 (t)] +M 2 (t) 
9b L J 

M 4 (t) = M 4 (t)+M 3 (t) + Mi(t) 



(39) 



The total magnetization is 



M tot {t) = J2 M i(t) = 2 



9B + 9A 
9B 



Mi{t) 



where the complex susceptibility x(w,T) is given by 



, 9b ~9A 

9B 



M 3 (t)=x(u)H t 



(40) 



a* 2 r 1 

x(w,T) = N — 1 - 71 771 - 72^2 + 2727i(*7i + %) 



.9b - 9 A \ 2 1 1 



A4 1 — 



(^) 



A3 1 - ILUT 3 



(41) 



The approximate static susceptibility, calculated from Eq. (|4"Tj) for zero frequency, is shown in Fig. 4a for I < r < 1. 
One immediately notices that, in striking contrast with the case < r < \ displayed in Fig. 3, the low temperature 
behavior of the static susceptibility is correctly reproduced. 

The latter feature appears at odds with the expectation that a decoupling approximation should work better the 
higher the temperature. However it is worth noticing that, for the ID ANNNI model, the T — > asymptotic behavior 
of the static two-spin correlation functions is very different depending on the value of r. For < r < 5 both the 
inter- and the intra-sublattice spin-spin correlations are strong (771 sa 772 ~ V3 ~ 1> see Note|32| later). In contrast, for 
r > i the intersublattice correlations are strong (772 ~ —1, sec Eq. (|45|) later), whereas the intrasublattice correlations 
are exponentially vanishing (77! w 773 w 0, see Eq. (1451) ). At intermediate temperatures intrasublattice correlations 
become significant, too, and the decoupling approximation becomes less satisfactory; at high temperatures, it works 
well again, since all correlations (both intra- and inter-sublattice) decrease. 

It should be remarked that the above considerations about the behavior of static correlation functions can not, on 
their own, account for the good agreement found, at low T, in the case h < r < 1. In fact, the use of equilibrium 
values for the spin correlations might be questionable, since the T — Glauber dynamics does not lead to the ground 
state of the ID ANNNI model in the entire region < r < l. 25 To this regard, first we observe that the physical 
mechanism which at T = prevents the system from reaching the ground state is different, for \ < r < 1, with 



respect to the case < r < 



1 25.27.28 



Next, considering that at T = a ID model is simultaneously in the ordered 



phase and at its critical point, while our theory applies at T > 0, we believe that some insight into the problem might 
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kT/J, 



FIG. 5: (color online) Temperature dependence of the amplitude of the complex susceptibility |x(a;,T)| for an ANNNI chain 
with Ji — 1, J2 = —1.25 (r = 1.25), subject to a weak external magnetic field oscillating at frequency u>. Figures (a), (b) refer 
to the compensated case (gA = ffs = 2), while Figures (c) and (d) to the uncompensated case (gA — 2, qb — 3), for selected 
values of the frequency (u/a — 0.001 and 10). In Figs. (c,d) the thin (color) lines represent the contributions to the amplitude 
of the two terms on the r.h.s. of Eq. (|41|l . while the thick (black) line is their sum. No resonant behavior is observed. 



be provided by a careful study of the role of a small but non-zero temperature on the coarsening of the ID ANNNI 
modelM 

In Fig. 4b the approximate static susceptibility, calculated from Eq. (|4"Tj) for zero frequency in the case r > 1, is 
reported. A nice overall agreement with the exact transfer matrix resul t 21 ' 22 is obtained. In this case our approximate 
results are expected to be quite reliable since the long-time approximation is well founded (for r > 1, the static 
equilibrium state can asymptotically be reached even at T = 0^ and thus the use of static spin-spin correlation 
functions is justified); moreover, the decoupling approximation is expected to be satisfactory both at high and low 
temperatures. Finally it is worth mentioning that, in the limiting case 1/r = (i.e., J\ = 0), the transfer matrix 
result for the static susceptibility is exactly reproduced by Eq. (j4"T]) for u> — (not shown). 

In Fig. 5 the temperature dependence of the amplitude of the complex susceptibility |x(w,T)|, obtained from 
Eq. (|41[) , of an ANNNI chain with nnn antiferromagnetic coupling dominating over the nn ferromagnetic one (com- 
petition ratio r = 1.25) is reported - for selected values of the oscillation frequency ui of the external magnetic field 
- both in the compensated (gA = Qb = 2) and uncompensated (gA ^ 9b) case. No resonant behavior was observed 
even in the uncompensated case since, in the T — ► limit, both the zero-field static susceptibility and the relaxation 
times (T3 and T4 in Eq. (T4"l"j) ) fail to diverge. Thus, for low frequency, a resonance condition - similar to the one in 
Eq. ([2U)l - can not be fulfilled. In the case \ < r < 1 a qualitatively similar behavior for |x(w,T)| was found (not 
shown) . 



C. Critical dynamics of the ID ANNNI model for r > 1 



The identification of r = 0, r = 1 and 1/r = as dynamic critical transition points for the ID ANNNI model 
with single spin- flip Glauber dynamics was recently proposed in theoretical studies of T = coarsening 2 ^- (i.e., the 
relaxation of the system into the ground state after a quench from high temperature) and T = persistence 2 ^ (i.e., 
the probability for a spin to remain in its original state after a quench from high temperature). In such T = studies, 
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the dynamic critical exponent z is customarily denned as the inverse of the growth exponent n of the domain size 

L(t)~t n -t 1 '*' . (42) 

For the ID nn Ising model, analytical calculations 2 ^ provided z' = 2. For the ID ANNNI model with r > 1 numerical 
calculation s 27 ! 30 predicted a somewhat higher dynamic exponent, z' ~ 2.3. Finally, it is worth noting that Sen and 
Dasgupta^ in their study of t = persistence in the ANNNI chain, found that the dynamic critical exponent z' 
undergoes abrupt changes for r — (when a slight amount of nnn interaction is added to the nn one), for 1/r = 
(when a slight amount of nn interaction is added to the nnn one), as well as for r = l£± 

The fair accuracy of our approximate theoretical approach in describing the low temperature static susceptibility 
of the ANNNI chain with r > 1, see Fig. 4b, encouraged us to tentatively estimate the dynamic critical exponent. 
However, since we work at finite temperature, rather than at T = 0, we use a different definition, namel y 11 ' 31 

-L-Wir (43) 



an V£ 

where Ai is the smallest eigenvalue of the dynamical matrix, see Eq. (f36|) . and £ is the static correlation length of the 
infinite system (the lattice constant c along the chain was set to 1). For the compensated case gA = SB, the latter 
quantity can be analytically calculated using the transfer matrix method^ see Eq. (IA18|i . and for r > 1 its expansion 
in the T — > limit turns out to be 

/1\2 1 2(J 1 -2|J 2 |) 

(-) *-e ** , (44) 

where we have explicitly taken into account that J\ > and J2 < 0. 

Taking into account the T — ► asymptotic behavior, for r > 1/2^ of the 7$ and 77, coefficients 

„ „ 4|J 2 

71 f» 1 — 2e *r 72 » — 1 + 2e 

1 Jj-g] JgJ £l— 2JJ3J 3 Jl-2|J 2 I , , 

»7i ~ 2 e kT ?72~-l + e fcr 773 w --e fcT , (45) 

for the inverse of the longest relaxation time we obtain, provided that J\ 7^ 

1 J1-2U2I , A . 

= Ai « e ^ . (46) 

an 

In the special case Ji = (i.e., 1/r = 0), letting 7! = in Eq. (j3"o| and using the T — > expansion for 72 in Eq. (|45[) . 
we obtain 

=Ai«2e— s#\ (47) 

In conclusion, within our approximate theoretical scheme, the dynamic critical exponent of the ID ANNNI chain 
with competing nn and nnn exchange interactions was found to be z — 1 for any finite r > 1, while in the absence 
of competing interactions (i.e., for r = and 1/r = 0) we found z = 2. Notice that, for the ID Ising model with 
exchange limited to the nn (r = 0), the value z = 2, obtained using the definition in Eq. ((43|) ^> coincides with the 
value z' = 2, obtained using the definition in Eq. (|42|).— This appears not to be the case for the ID ANNNI model 
with 1 < r < +00, where the values z = 1 (present work) and z' ~ 2.3 (References [27l|30t) were found. In order to 
ascertain the origin of this discrepancy, we believe that it would be useful to study the role of a small but non-zero 
temperature (T > 0) on the coarsening dynamics of the ID ANNNI models 



IV. CONCLUSIONS 



In conclusion, in this paper we have studied the effect of antiferromagnetic interactions on the single spin-flip Glauber 
dynamics of two different one-dimensional (ID) Ising models with spin ±1. For the first model, an Ising chain with 
antiferromagnetic exchange interaction limited to nearest neighbors and subject to an oscillating magnetic field, the 
system of master equations describing the time evolution of sublattice magnetizations can easily be solved within a 
linear field approximation and a long time limit. Resonant behavior of the magnetization as a function of temperature 
(stochastic resonance) is found, at low frequency, only when spins on opposite sublattices are uncompensated owing 
to different gyromagnetic factors (i.e., in the presence of a ferrimagnetic short range order). For the second model, 
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the axial next-nearest neighbor Ising (ANNNI) chain, where the nnn antiferromagnetic exchange coupling is assumed 
to compete with the nn ferromagnetic one, the long time response of the model to a weak, oscillating magnetic field 
is investigated in the framework of a decoupling approximation for three-spin correlation functions, which is required 
to close the system of master equations. Within such approximate theoretical scheme, the T = dynamics of the 
Ising-Glauber chain with competing interactions is found to be in a different universality class than that of the Ising 
chain with antiferromagnetic exchange limited to nearest neighbors (r = 0) or limited to next-nearest neighbors 
(1/r = 0). In particular, we find an abrupt change in the T = dynamic behavior of the model in the neighborhood 
of the dynamic critical point 1/r = since, when a slight amount of ferromagnetic nn exchange is added to the 
antiferromagnetic nnn exchange, we find that the critical exponent z, defined by Eq. (|43|) . changes abruptly from 
z = 2 to z = 1. Considering that z — 2 is also the value of the dynamic critical exponent for the unfrustrated nn Ising 
chain, one might expect similar abrupt changes in z to occur also in the neighborhood of the dynamic critical points 
r = (i.e. when a slight amount of AF nnn exchange is added to the nn F exchange) and r = 1, as suggested by 
studies of T = coarsening dynamics^ and T — persistence^! in the ANNNI chain. Unfortunately, the inaccuracy 
of our approximate theoretical scheme in reproducing the static susceptibility of the ID ANNNI model with < r < 1 
for low temperature prevented us from calculating the dynamic critical exponent in this range of the competition 
ratio. 



APPENDIX A: ANALYTIC TRANSFER MATRIX RESULTS FOR THE STATIC PROPERTIES OF ID 

ISING MODELS 



1. The ID nearest neighbor Ising model with alternating spins in a static field 



In this subsection we calculate, within the transfer matrix formalism^ the static properties of the ID Ising model, 
Eq. (Q}, with nearest neighbor coupling J of either sign, subject to a static magnetic field H (i.e., u — 0). Two types 
of spins with different gyromagnetic factors (gA ^ 9b) are assumed to alternate along the chain. Taking periodic 
boundary conditions, the partition function of the chain of length TV (with N even without loss of generality) can be 
expressed as 



Zn = Tr(e St ) = ^ " ' K (ai, a 2 )L(a 2 , a 3 ) ■ ■ ■ K(a N ^i, cr N )L(a N , <ti) 

(7l=±l (T2— ±1 (TiV— ±1 



(Al) 



where, letting J = J/(kT), Ha — (g aHqHq) / '(kT) , h B = (gB^oHo)/(kT), the two different kernels K and L are 
defined as 

K(a 2i -i,(T2i) = e Ja2i -^ i e^ hA ' T2i - 1+hB ' J2l) L(a 2i , <r 2 i+i) = e Ja2ia2i + 1 e^ hB<T2i+hAa * i + x \ (A2) 
Summing over the even sites, Zjq can be expressed as 

W JV 

Z N = (A+) 



(A- 



(A3) 



in terms of the eigenvalues 

A± = e 2J cosh(fiA + h B ) + e~ 2J cosh(hA - h B )± 



± J e AJ cosh 2 + h B ) + e~ AJ coslr (/^ - h B ) + 2 cosh(/i A + h B ) cosh(h A - h B ) + 2 



e 4 ^ - e-^ (A4) 



of the real symmetric 2x2 matrix 

a 2J+h A +h B 

S= ' 



^2J-h A -h[ 



-2J-h A +h B 



(A5) 



It is immediate to verify that, in the limit gA = g B , the well-known result for the ID nn Ising chain in a static external 
field is recovered^ In the thermodynamic limit N — > oo, only the larger eigenvalue A + matters, Zn — + (A + )~, and 
the static susceptibility in zero field can be expressed in terms of its second derivative with respect to the field H 



x(w = 0,T) 



N , 



-kT 



1 d 2 A, 



A+ dH 2 



H=0 



kT 



e kT + 



9 a - gB 



e kT 



(A6) 
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2. The ID ANNNI model in zero field 



In this subsection we collect, for the reader's convenience, some exact results for the static properties of the ID 
ANNNI model in zero field, Eq. (|2"2")l , which were obtained by Stephenson-! and Harada22 in the case of a linear chain 
with N identical spins (gA = 9b = 9 an( l a = Using the transfer matrix method, the partition function can be 
exactly expressed as 



Z N = (A+) iV + (A_) 
in terms of the eigenvalues of the symmetric 2x2 matrix 



S = 



a c 
c b 



■'2 + -'l 



N 



g kT £ kT 
e kT e kT 



The eigenvalues take the formal 

A± = 

where 



1 -'2 

-[a + b±A] =e~^T 



4J 2 



A = v/(a-6) 2 + 4c 2 = 2e^ 



(A7) 



(A8) 



(A9) 



(A10) 



Both A+ and A are always real positive quantities. 

In the thermodynamic limit N — > oo, the static two spin correlation function rj n take the formal 



1 



Vn 



2 (A + r 



1 



AA' 



1 - 



a 2 -b 2 
AA' 



where the quantities A', defined as 



A' = y/( a + bf - 4c 2 = 2e^ J cosh 2 -e"W, 



and 



j£2. 



sinh 



■(£) 



( J\ \ 4J 2 
I I — C kT 



± \ / cosh' 



\kT) 



(All) 



(A12) 



(A13) 



may be complex. More precisely, the quantities \x± are real for T < To and complex conjugates for T > Tjj. Tjj is 
the so-called disorder point, defined by the equation A'(Td) = 0, which has solutions for < r < 1/2 at some finite 
temperature Tjj. For T < Tjj the static equilibrium two-spin correlation functions rj n = {(jj<jj +n ) present a monotonic 
exponential decay, while for T > Tu they have an oscillating exponential decay^ 

Summing over all pair correlations, the exact zero field static susceptibility can be expressed as2i 



The wave-vector dependent susceptibility, defined as 



a(a - b + A) + 2c 2 



b{b - a + A) + 2c 2 



kT 

n 

presents a maximum at a wave- vector g m , which is given by^ 

(m+ + m-)(a+-a_) 



cos q„ 



(A14) 



(A15) 



(A16) 



For < r < 1/4 one has q m — at all temperatures , while for 1/4 < r < 1/2 there is a definite temperature Tj, 
(y£ Td) above which q m ^ 0, whereas for T < Tl q m = 0. When 1/2 < r, one has q m {T = 0) = 7r/2. In the limit of 
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T — > oo, <? m tends to the mean field value cosq m = l/(4r). Expanding x(l) in the neighborhood of q m up to second 
order in Aq — q rn — q, one obtains a Lorentzian form, and the correlation length £ can be defined in terms of its full 
width at half maximum as 

= i x i g rL 2 ( ai7 ) 



and turns out to be 



i + e 2 (A g ) 2 



fl\ 2 _ (A+ - A_ - n+ - n-f 
\0 - (m + +/x-)(A + -A_)-4 M+ M- t0r<?m - U ' 
1x2 ( Ai+ - M _) 2 (A+ + A_) 2 



^ ( M+ + Ai _) 2 (A + -A_) 2 -16 / i 2 / i 2 _ 



for q m + 0. (A18) 



For 1/4 < r < 1/2, it turns out that at Tl the correlation length becomes zero, which is a characteristic of the Lifshitz 
points 
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